Insulating state of granular superconductors in a strong-coupling regime 
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We analyze the possibility of the formation of a magnetic-field-induced insulating state in a two- 
dimensional granular superconductor with relatively strong intergranular coupling and show that 
such a state appears in a model with spatial variations of the single-grain critical magnetic field. 
This model describes realistic granular samples with the dispersion in grain sizes and explains the 
mechanism leading to a giant peak in the magnetoresistance. 
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I. INTRODUCTION 

Recent experiments on superconducting films revealed 
a giant peak in the low temperature magnetoresistance 
with the growth from several times (in polycrystalline 
samples 1 ) to several orders of magnitude (in amor- 
phous films 2 ^^), as compared to the value near the 
superconductor-insulator transition, and the re-entrant 
drop upon further increase of magnetic field. The 
magnetic field-induced insulating phase and the related 
superconductor-insulator transition are long known, and 
extensive theoretical and experimental effortsSi^ re- 
sulted in a remarkable progress in understanding many 
features of these phenomena. Yet the observed giant peak 
in magnetoresistance is puzzling and poses a challenge for 
theory. Theoretical modelsi2*ii suggested possible devel- 
opment of spatial inhomogeneities in amorphous samples 
and indicated the role they may play in formation of the 
resistivity maximum. 11 Consequences of the existence of 
the superconductor islands network were explored by the 
numerical studies in Ref. IT^ : the computed behavior of 
magnetoresistance was in a qualitative agreement with 
the experiment. 

In the present paper we propose and investigate a 
model for polycrystalline and/or granular films which of- 
fers a mechanism leading to the giant peak in the mag- 
netoresistance. We model the superconducting film in 
a strong magnetic field as an array of alternating super- 
conducting (S) and normal (N) granules and demonstrate 
that the superconducting insulator (SI) phase appears in 
the strong coupling regime. In real multi-disperse films 
the formation of a network of S and N grains occurs nat- 
urally as a consequence of the spatial variations of the 
single grain critical fields that inevitably appear due to 
dispersion of grains' sizes. Magnetic field close to the 
average single-grain critical field drives about a half of 
the grains into a normal state, giving rise to the network 
of S and N granules, see Fig. ^ Even at a strong inter- 
granular coupling, this array is subject to the Coulomb 
blockade (similar to a single SN junction^), leading to 
the insulating state of the sample. Varying the magnetic 
field results in the change of the relative fraction of the 
S and N grains. 

Note that the "standard" model of an array of iden- 
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FIG. 1: (a) Superconductor film consisting of grains of slightly 
different sizes, (b) The same granular array under applied 
magnetic field close to the average single-grain critical field. 
Varying magnetic field results in a change of the relative 
concentration of superconducting (white) and normal (gray) 
grains. 

tical S granules^ cannot describe the formation of an 
insulating state in the regime of a strong intergranu- 
lar coupling. The charging energy of a single grain is 
renormalized 14 down to the values E c ~ A/G where A is 
the order parameter in a single grain and G is the inter- 
granular tunneling conductance. The Josephson coupling 
Ej ~ GA ^> E c if G 1, and this results in the super- 
conducting rather than the insulating stated 

We consider the case where most of the grains are so 
small that the Zeeman effect is at least of the same or- 
der as the orbital one, and, as a result, the transition 
between the superconducting and normal states within 
a single grain is of the first order>i£ Such samples are 
the good candidates for demonstrating the SI phase since 
the applied strong magnetic field H close to the Clogston 
field H will form a strongly inhomogeneous mixture of S 
and N grains with the order parameter in the S granules 
close to that in the absence of the magnetic field, Ao. 
Figure [3 shows the phase diagram of a single supercon- 
ductor grain (the derivation is presented below) drawn 
in the coordinates Zeeman energy h = [IbH, where hb 
is the Bohr magneton, vs. the size-dependent orbital 
effect pairbreaking parameter a [in spherical grains 1 ^ 
a = D(eHa/2c) 2 /5, where a is the grain diameter and D 
is the bulk diffusion constant] . Varying the grain size one 
crosses the line of the phase transition which at strong 
enough fields is indeed always of the first order. 

We now recall a well known result of the percolation 
theory that in the site problem for the two-dimensional 
black and white alternate square (or honeycomb) lattice 
there exists a range of the relative concentrations of the 
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FIG. 2: Phase diagram of a superconducting grain at T = 
in coordinates Zeeman energy h vs. pairbreaking parameter 
a. The dashed line separates the gapless (GS) and gapful (S) 
superconducting states. First- and second-order phase transi- 
tion lines are shown by thick and thin solid lines, respectively. 



sites where simultaneously neither black nor white sites 
percolated In our model, such a "no-percolation" state 
inevitably leads to transport via sequences of SN junc- 
tions; this is an essential ingredient for explaining the SI 
phase. To capture this qualitative fact and at the same 
time to simplify the calculations, we consider a system 
with the checkerboard arrangement of the S and N gran- 
ules. 



A. Main results 

We show that the applied strong magnetic field, close 
to the Clogston field H — A / v2/i b , causes the forma- 
tion of an insulating state with the gap in the electron 
excitation spectrum given by 
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where G is the tunneling conductance^ (in units of 
2e 2 /h) and E c is the single grain charging energy. This 
suggests the conductivity behavior 
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Equation is obtained in the approximation of the 
nearest neighbor tunneling between the grains. 

To verify whether the higher-order tunneling processes 
can indeed be neglected, we analyze (i) electron tunneling 
between two N grains via a virtual state in an S grain 
and (ii) Cooper pair tunneling between S grains via an 
N grain. Small amplitudes of these processes ensure the 
stability of the insulating state with respect to formation 
of the metallic and superconducting state, respectively. 
Analyzing the corresponding corrections we show that 
the insulating state is stable as long as G < G* , where 
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where <5 is the mean energy level spacing in a single grain. 
Since the ratio A /<5 is large (otherwise even the mean 
field approximation for a single grain cannot be used), 
G* is larger than the critical value 



r ~ 1 i Ec 
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that marks the onset of strong Coulomb correlations in 
the two-dimensional granular array in the absence of su- 
perconductivity^2*2S Thus, for samples with intergranu- 
lar conductance laying within the interval 



G c < G < G* 



(5) 



the insulating state is indeed induced by the local super- 
conducting correlations. 



II. PHASE DIAGRAM OF A 
SUPERCONDUCTING GRAIN 

Now we turn to the analytical derivation of our results. 
We begin with the discussion of a single granule of the 
size less than the coherence length. The difference be- 
tween its thermodynamic potentials in superconducting 
and normal states can be written as 
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where A is the electron-phonon coupling constant, Co = 
to — ih, lu is the fermionic Matsubara frequency, and g, /, 
and p are the Usadel Green functions^ subject to the 
constraint g 2 + ff^ = 1. Varying f2 with respect to f> un- 
der the above constraint results in the zero-dimensional 
Usadel equation 



(7) 



while varying in A* yields the self-consistency equation 
A = XnTj2f, (8) 



(below we choose A real and / = f'). The second-order 
phase transition line can be obtained by solving Eq. (JTJ 
in the first order with respect to A and inserting the 
resulting / function into Eq. ©. At zero temperature, 
T = 0, we arrive at2S 
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The second-order phase transition line turns into the 
first-order one (see Fig. 0) that can be found by solving 
Eqs. and © under the condition fi = 0. At a = 
we reproduce the Clogston result 15 for the critical field 
ho = Aq/x/2. 
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The boundary between the gapful (S) and gapless (GS) 
phases is obtained from the condition that the gap in the 
electron spectrum A„ = QiML2£ 
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where A has to be found from Eq. JSJl. The latter is not 
affected by the Zeeman term (since the spin susceptibility 
in this regime is zero at T = 0) in the gapful regime and 
thus assumes a simple form^ 



A ira 
ln X = 4A- 
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We see that at small enough grain sizes (small a) the 
gapless region does not exist. In what follows we will 
focus on the situation where the superconducting grains 
all have the gap. The validity of this assumption is sup- 
ported by the fact that superconducting insulator phase 
was observed at very high magnetic fields of the order of 
the Clogston limit. The gapless region is not favorable 
for the formation of the SI state. 



III. ARRAY OF SUPERCONDUCTING GRAINS 

We describe the array of S and N grains by the phase 
actionSi 



(12a) 



where the first term is the Coulomb part of action with 
the sum going over all grains. The term S ns describes 
the coupling of normal and superconducting grains: 
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where summation goes over the nearest neighboring 
grains, the phase difference 4>ij = 4>ij{r) — </>i(r) — (f>j(r) 
and the kernel 

a(r) =T 2 Y e-^-^ T g s {u)g n {u') (13) 



is expressed through the Usadel Green functions g s (u) 
and g n (to) of the S and N grains. 

In the limit G>G C the term S ns can be simplified by 
expanding it up to the second order in </>y (t) and using 
the local approximation for the kernel a to 
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where g' s (0) is the derivative of the function g with re- 
spect to u> taken at to = 0. It can be found from Eq. J7J, 



g' s (0) = (-a + A// 3 (0)) 1 , and in the limiting case that 
we assume, h w h ^ a, reduces to 
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In the Fourier representation we obtain 
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where we have introduced the notations 
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with a being the lattice vectors, and 
B = SGg'M- 



(17) 
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Integration over quasimomentum q is over the first Bril- 
louin zone. Equation (|16|) can be interpreted as the renor- 
malization of the charging energy 
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The gap in the electron spectrum becomes 
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where Eqq is a diagonal element of the Coulomb inter- 
action matrix that is given by the inverse Fourier 
transform of E c (q). Within the logarithmic accuracy 
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leading to Eq. in the case of strong fields where Eq. 
(|T5|l can be used. We see that A^ e ^ is indeed less than 
Ao as long as G ^ G c . This justifies the local in time 
approximation used in obtaining Eq. (|14|l . We note that 
in terms of the phase functional there is no difference be- 
tween the S and N grains. However, one has to remember 
that the S grains have the "big" single particle gap A g . 
At small nonzero temperatures, the electron excitations 
"live" only in the N grains, above the gap . Despite 
the superconducting gap, electrons can tunnel between 
two N grains via a virtual state of an S grain, see dia- 
gram in Fig. |HIa) (this higher-order tunneling process is 
analyzed below) . Therefore the transport has the activa- 
tion form with the gap A^ e ). 

Now we consider the higher-order tunneling processes 
so far neglected and find the conditions, at which the 
contribution of the diagrams in Figgis not essential. 



A. Stability of the superconducting insulator state 
with respect to the formation of the normal state 

The diagram in Fig. 01a) should be analyzed in or- 
der to check the stability of the SI state with respect to 
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FIG. 3: Diagram (a) represents the coupling of the normal 
grains (Ni and N2) via the superconducting grain S. Dia- 
gram (b) describes the processes of Cooper pair tunneling 
between the superconducting grains (Si and S2) via the nor- 
mal grain N. Solid (electron) lines are connected via dashed 
lines that represent impurities, forming the zero-dimensional 
diffuson/Cooperon. Solid lines with two arrows represent the 
anomalous Green functions of the superconductors. 



B. Stability of the superconducting insulator state 
with respect to the formation of the superconductor 

state 

Now we analyze the stability of the insulating state 
with respect to formation of superconductivity that may 
be established due to the Josephson coupling of the S 
grains, see the diagram in Fig.^b). To this end we in- 
troduce an effective Hamiltonian acting on Cooper pairs 
and defined on the sublattice of the S grains 



H 



J 



(28) 



formation of the normal state. Such tunneling processes 
lead to an effective coupling of N grains that due to con- 
dition A^ < A can be viewed as an effective tunneling 
coupling with the conductance 2 ^ 



(22) 



2 and 



where the dimensionless matrix element tij 
tij = 1 for nearest and next-to-nearest neighbors in the 
N sublattice, respectively. This effective coupling results 
in the additional term in the phase action 
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To find the correction to the Coulomb gap A^ e ) in the 
first order in S nn , we can consider the normal sublattice 
of the array only. The problem then becomes equivalent 
to the calculation of the correction to the Mott gap in the 
array of N grains. Using the result of Ref. |2(] we obtain 
the correction to A^ e ^: 
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Using that E C G > A from Eqs. lfT§|> and (|2~1J) we obtain 



(25) 



Comparing Eqs. and (|25|) . we conclude that <5A( e - ) <C 
A( e > if 
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We note that the effective coupling IL'l't also defines the 
prefactor in the conductivity behavior: 



4e 2 / A< e ) ' 

a=^G 2 ^(0)exp 



where the function g' s (0) was defined in Eq. I|15|l . 



(27) 



The Hamiltonian is written in terms of the superconduct- 
ing phase tp, which is related to the "normal" phase <f> as 
<p = 2cf). Here i, j belong to the S sublattice, n = —id/ dip 
is the Cooper pair density operator, J is the effective 
Josephson coupling of S grains via a N grain, and tij — 2 
and tij = 1 for nearest and next-to-nearest neighbors 
in the S sublattice, respectively. The insulating state of 
the model l|28ll is characterized by the energy gap in the 
spectrum of Cooper pair excitations^ Due to the dou- 
bled Cooper pair charge in the limit J = this gap is 



A (s) = 4A (e) . 



(29) 



The Josephson coupling J leads to the suppression of 
that in the limit of weak coupling can be found via 
the straightforward perturbation theory 26 



<5A« = -^^t M = -6J. 



(30) 



The condition SA^ < A^ is satisfied if J < A^ s ) ~ 
A / G. The diagram in Fig. |3^b) in the case of the large 
Zeeman energy, h w ho 3> a, gives^i 



(31) 



and we arrive at the estimate for the upper boundary 
for the conductance One can verify that the two- 

electron (Andreev) tunneling processes between the nor- 
mal grains via an intermediate superconducting one are 
also suppressed under the same condition. 



IV. COMPARISON WITH EXPERIMENT 

Now we discuss the relation between our results and 
experimental data of Refs. Q and 0- A note in order is 
that the proposed model of a granular superconductor 
in the vicinity of a superconductor-insulator transition 
as an normal/superconductor granules alternating chess- 
board seems to capture adequately general features of 
the systems of real morphology and, therefore, allows for 
the comparison with the experimental data within its 
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range of applicability. Our results hold for conductances 
G > G c , where the critical conductance G c is given by 
Eq. (@J. This implies that the value of a film conduc- 
tance G may be less than unity in real experiment. The 
resistance of the samples used in experiments of Refs. Q] 
and0 was of the order of 20 kfi, which corresponds to 
the dimensionless conductance G ~ 0.1, i.e., falling on 
the boundary of applicability of our theory. Yet one can 
hope that our results can be used as a crude estimate. 

In Ref. Q], the data showing a noticeable peak in mag- 
netoresistance were presented only for the high resistivity 
samples. The maximum resistance does not demonstrate 
the activation dependence on temperature. Apparently 
our theory does not apply to this case, presumably due 
to low intergranular tunneling conductance. 

Next, our theory is developed for granular systems and 
therefore should not apply to homogeneously amorphous 
samples. However, it is possible that realistic samples 
do not exactly correspond to one of the two limiting 
cases. For example, amorphous samples can be not uni- 
formly disordered but rather consisting of separate re- 
gions (due to fluctuations of chemical composition or im- 
purity concentration) favoring superconductivity. In this 
case the system reminds the granular samples studied in 
the present paper. 

Keeping this possibility in mind, we make estimates 
for the experiment of Ref. 0, where amorphous InO sam- 
ples were investigated. We note that in the case of InO 
films the formation of the "superconducting insulator" 
state was observed in the samples where the values of 
the resistance at high magnetic fields, while falling into 
the metallic domain, were close to the critical resistance 
that separates insulating and metallic phases at T — > 0. 
Thus, the estimate of the gap A^ e ) in Eq. JIJ requires 
taking G » G c . Plugging this into Eq. QJ, we obtain the 
gap 
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Within the logarithmic accuracy, the ratio of the loga- 
rithms is about unity. Finally, in order to compare our 
estimate with the experimental findings, we express the 
superconducting gap A in the absence of the magnetic 
field through the critical temperature T c . Using the stan- 
dard relation Aq = 1.76 T c , we arrive at the estimate 
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or, equivalently, T c j A*- e ) ss 3, which certainly is in a rea- 
sonable agreement with the experiment. Indeed, accord- 
ing to Fig. 4 of Ref. 0, the ratio T c /A^ depending on 



the sample varies between 0.6 and 3. Bearing in mind 
that our model was strictly speaking designed for granu- 
lar samples, one can hardly expect a better agreement. 
V. DISCUSSION AND CONCLUSIONS 

Here we comment on the conductivity behavior in the 
vicinity of the magnetization peak: At lower magnetic 
fields the concentration of the superconducting compo- 
nent increases resulting in the formation of the direct 
Josephson coupling between the superconducting grains. 
The Cooper pair gap decreases and eventually be- 
comes smaller than the electron gap A^ e '. At the same 
time there is no reason for the strong suppression of the 
electron gap due to the Josephson coupling. Thus, closer 
to the insulator to superconductor (IS) transition at finite 
but low temperatures, the transport is mediated by the 
activation of Cooper pairs rather than by the single elec- 
tron excitations, hence a ~ exp (-A^ /T) at T < A( s ). 
Thus, we expect that the resistivity within the IS transi- 
tion may be described within the effective Bose model as 
in Ref. 0. At the opposite side of the magnetoresistance 
peak (i.e., at higher fields) the fraction of the normal com- 
ponent increases leading to the direct coupling of normal 
grains. This results^ in the decrease of the electron gap 
A^ 6 ) and eventually in the destruction of the insulating 
state. Note that for the SI state which we studied the 
spatial orientation of the applied magnetic field was not 
essential, however closer to the IS transition one expects 
it to become important. 

In conclusion, we have investigated formation of an in- 
sulating state (corresponding to a giant peak in the mag- 
netoresistance) in a two-dimensional granular array with 
alternating superconducting and normal granules. We 
have demonstrated that this model is a good represen- 
tation for realistic granular (or polycrystalline) systems 
with the grain size dispersion in a strong magnetic field. 
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